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Clustering at high redshift: The connection between 
Lyman Alpha emitters and Lyman break galaxies 
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ABSTRACT 

We present a physically motivated semi- analytic model to understand the clustering of 
high redshift Lyman Alpha Emitters (LAEs). We show that the model parameters con- 
strained by the observed luminosity functions, can be used to predict large scale bias 
and angular correlation function of LAEs. These predictions are shown to reproduce 
the observations remarkably well. We find that average masses of dark matter halos 
hosting LAEs brighter than threshold narrow band magnitude ~ 25 arc ^ IO^^Mq. 
These are smaller than that of typical Lyman Break Galaxies (LBGs) brighter than 
similar threshold continuum magnitude by a factor '^ 10. This results in a smaller clus- 
tering strength of LAEs compared to LBGs. However, using the observed relationship 
between UV continuum and Lyman-a luminosity of LAEs, we show that both LAEs 
and LBGs belong to the same parent galaxy population with narrow band techniques 
having greater efficiency in picking up galaxies with low UV luminosity. We also show 
that the lack of evidence for the presence of the one halo term in the observed LAE 
angular correlation functions can be attributed to sub-Poisson distribution of LAEs 
in dark matter halos as a result of their small halo occupations. 

Key words: cosmology: theory - cosmology: large-scale structure of universe - galax- 
ies: formation - galaxies: high-redshift - galaxies: luminosity function - galaxies: statis- 
tics - galaxy: haloes 



1 INTRODUCTION 

Over the past decade there has been a growing wealth of 
observations probing the properties of high redshift galaxies. 
Various surveys, using th e Lyman break color selection 
technique (see for example. iMadau et al.ll 19961: ISteidel et al.l 



Lyman-a emitters (LAEs). These surveys helped to infer 
the statistical properties of LAEs, particularly their UV 
and Lyman -a luminosit y funct i ons an d angular corr e lation 
functions jOuchi et al] |2003|. l2008l : IShiova et all [200i; 



"'^P ■'.': ^r~ 'i ~ i 

Il996l : lAdelberger et aLlll998l : ISteidel et al.lll998l '). provided 



IWang et al 



20091 : lOuchi et al.lbOlOl ) While LBGs selected 



fairly g ood estimates of luminosity functions (LF) u p to 
z ~ 8 (|Bouwenset al.ll200l l2008l : iReddv et al1l2008l) and 
spati al clustering u p to z ~ 5 (IGiavalisco fc DickinsonI 



20011: lOuchi et al.l l200i. l2005l: [Kashikawa et all l2006l: 



Hildebrandt et all l2009l : ISavov et all l201ll : iBielbv et all 
20 111 ) of these Lyman break galaxies (LBG). On the other 
hand, narrow-band s ea rches for hig h reds hi ft Lym an-a 



line (iHu et al 

20031: Shimasaku et al 



19981: iRhoads et al.i I2OO0I: 



200€- 



Gronwall et al 



Ouchi et al.l 
20071: 



Finkclst ein et al.l 1 2007: D awson et a l. 2007; I Ouchi et all 
2008: Shiovaet al.li2009 : Wang et al.i,2009 : Hu et al.ll2oTa 
Ouchi et al.l I2OI0I ) detected have substantial number of 



through colour cuts are biased towards bright galaxies 
the narrow band selection of LAEs is biased towards 
galaxies having strong emission line and weak continuum 
emission. Therefore, these two techniques pickup galaxies 
with different type of selection biases. Despite several 
detailed studies there is no clear consensus on whether 
there is any differences between the two populations based, 
on properties such as stellar mass, dust co ntent, age and 
star formation rate etc (see for example , iGawiser et al 



20061: iPentericci et al.ll2007l : lOuchi et alll200a : iKornei et al 



20101 ). 
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In the hierarchical model of structure formation the 
statistical properties of galaxies are determined by that of 
the parent dark matter halo population, given a prescrip- 
tion for how stars form inside these halos. The properties of 
dark matter halos are quite well understood using N-body 
simulations and analytical models like the halo model of 
large scale structure. There have been a number of stud- 
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ies over the past years, probing the properties of LAEs 
using semi-analytical a n d numerica l (simulat ion) methods 
llLe Delliou et al.1 l2006l: iMao et al.1 p 007; Ko bavashi et al.1 
20071: ISamuiet al booj lOaval et amOQg.: .Tilvi et al.ll20o'9F 



Samui et al.ll2009l: Zheng et al.ll201ol: iNagamine et al.ll2oTc 



Kobavashi et al.l2010l : lGarel etal.ll2012l '). These models have 
been successful in reproducing the (i) UV LF of LBGs, (ii) 
UV LF of LAEs and (iu) Lyman-a LF of LAEs. We also 
have been extensively developing simple physical models 
of galaxy formation to understand the LFs of LBGs and 
Lyman-a emitters. Any galaxy formation model that can 
simult aneously explain the ob served properties (like for ex- 
ample, iDaval fc Ferraral l2012h of LBGs and LAEs is very 
useful to further our understanding of physics of galaxy for- 
mation at high-z. 

In particular by modelling the clustering of LBGs and 
LAEs one will be able to understand any differences in 
the range of halo masses probed by these two population 
of sources. Recently we have presented a physically moti- 
vated semi-analytic halo model to simultaneou sly explain 
the L F and spatial clustering of high-z LBGs (|jose et al.l 
120131 ). Here we extend this model to explain the measured 
LFs and angular correlation function of high-z LAEs. The 
organization of this paper is as follows. In the next section 
we briefly outline our models for computation of the LF of 
LAE. In Section 3 we expand our models for LBG cluster- 
ing to explain the clustering of LAEs. Section 4 presents 
a comprehensive comparison of the total angular correla- 
tion function computed in various models with observa- 
tions. A discussion of our results and conclusions are pre- 
sented in the flnal section. For all calculations we adopt a 
flat ACDM universe with cosmological parameters consis- 
tent with 7 year Wilkins on Microwave Ani sotropy Probe 
(WMAP7) observations (JLarson et al.l 1201 J ). Accordingly 
we assume n„ = 0.27, ^a = 0.73, Q.b = 0.045, h = 0.71, 
Us = 0.963 and ag, — 0.801/i~^Mpc. Here Q,i is the back- 
ground density of any species 'i' in units of critical density 
Pa- The Hubble constant is Ho = IQQh km s^^ Mpc~^ 



2 THE STAR FORMATION RATE AND 
LUMINOSITY FUNCTION 

Here we recall the salient points of our semi-analytic model 
used to e xplain the observe d luminosity function of high- 
z LBGs JSamui et al.l [20071 . hereafter SSS07) and LAEs 
(|Samui et al.l |2009| . hereafter SSS09). The star formation 
rate (Msf) in a dark matter halo of mass M collapsed 
at redshift Zc and ob s erved at redshift z is given b y (see, 
IChiu fc Ostrikeill2000l : IChoudhurv fc Srianand|[20o3) , 



MsF^M, 



i Ij^ 



t{z)-t{z,) 

dyn(za) 



(1) 



X exp 



t{z)-t{Zc) 



where, the amount and duration of the star formation is 
determined by the values of /, and k respectively. Further, 
t{z) is the age of the universe; thus t{z) — t{zc) gives the age 
of the galaxy at z that has formed at an earlier epoch Zc, 
and tdyn is the dynamical time at that epoch. The stars are 
formed with a Salpeter initial mass function (IMF) in the 
mass range 1 — 100 Mq . The star formation rate of a galaxy, 



as given in Eq. ([T]), is then convolved to obtain the UV 
luminosity at 1500A (i.e Lisoo) and AB magnitude (A/as) 
at any given time using the procedure described in detail 
in SSS07. The observed luminosity is only a fraction l/r) 
of the actual luminosity because of the dust reddening in 
the galaxy. Having obtained the Mab of individual galaxies 
we can compute the luminosity function (^[Mab, z) at any 
redshift z using. 



^{MAB,z)dMAB 



^__ dn{M{MAB),Zc) dM 



dZc 



X -TTT — dMAB- 



dLir, 



dMAB 

Here dn{M,Zc)/dzc = h[M,Zc)dt/dzc, and h{M,Zc)dM is 
the formation rate of halos in the mass range (M, M + dM) 
at redshift Z c- We model this forma tion rate as the time 
derivative of ISheth fc TormenI (|l999l ) (hereafter ST) mass 
function as they are found to be good in reproducing the 
observed LF of high-z LBGs. Therefore we use h{M, Zc) — 
dnsT{M,Zc)/dt where nsT{M,Zc) is the ST mass function 
at Zc. Also note that we use the notation n{M) for dn/dM 
for convenience. 

Star formation in a given halo also depends on the cool- 
ing efficiency of the gas and various other feedback processes. 
We assume that gas in halos with virial temperatures {T-uir) 
in excess of 10* K can cool and collapse to form stars. In 
the ionized regions we incorporate radiative feedback by 
a complete suppression of galaxy formation in halos with 
circular velocity Uc ^ 35 km s~'^ and no su ppression with 
Vc^Vn = 95 km s"^ JBromm fc Loebll2002l ). For interme- 
diate circular velocities, a linear fit from 1 to i s adopted 
as the suppression fa ctor iBromm fc Loebl (|2002h : see also 



iBenson et all (|2002l ): iDiikstra et al.1 (|2004h . SSS07]. In ad- 
dition, we also incorporate the possible Active Galactic Nu- 
clei (AGN) feedback that suppresses star formation in the 
high mass halos, by multiplying the star formati on rate by 
a factor [1 -I- M/Magn\~°-^ (see I Jose et"alll2013l ). This de- 
creases the star formation activity in high mass halos above 
a characteristic mass scale Mgpn, which is be l ieved to be 
~ 10^^ Mq (see iBower et al.l 120061 : iBest et al.l l2006l). The 
crucial parameter of our model is ft/rj which governs the 
mass to light ratio of galaxies at any given redshift. This 
is fixed by fitting the observed UV LF of LBGs using x^ 
minimization. 

The Lyman-a luminosity produced in any star forming 
region (when we use the case-B recombination) is related to 
its instantaneous star formation rate and is given by, 



LLyc, = 0.67fe!/a(l " fesc)N-,MsF- 



(3) 



Here, hvc, = 10.2 eV is the energy of a Lyman-a photon and 
we fix fesc = 0.1, the escape fraction of UV ionizing photons. 
Further, A''.^ is the production rate of ionizing photon per 
unit solar mass of star formation which is a function of IMF 
and the metallicity. We use A^.^ = 10450 for the Salpeter 
IMF and metallicity oi Z = 0.0004 (See SSS07). 
The observed Lyman-a luminosity is given by, 



-' Lya 



rLya jr 
J esc J^Lya- 



(4) 



Here, fcJi°' is the escape fraction of the Lyman-a photons 
which is decided by the dust optical depth, velocity field of 
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Figure 1. Top panels show the predicted UV LF of LBGs (dashed hnes) at z = 3,4,6 and 7 a nd that of LAE s (thick Hne s) at 
z = 3.1, 3.8, 5 .7 and 6.6 along with the observed data . The observed data points of LBGs taken from lBouwens et al.l 1 120071. 120081) and 
[ReddvetaU l l2008r arc shown as red filled squares. The black filled triangles are the observed UV LF data points from lOuchi et alj 
1I2OIOI . 120081) . Bottom panels show our model predictions of Lyman a LF of LAEs at 2 = 3.1,3.7,5.7 and 6.6. Here, the observed data 



points and error bars are from lOuchi et alj ll20ld , I2OO. 



the ISM in the galaxies and the Lyman-a optical depth due 
to ambient intergalactic medium around the galaxies. 

Observations suggest that only a small fraction of LBGs 
(Gf) show detectable Lyman-a emission. In our models we 
assume that the UV LF of LAEs at any 2; is a fraction Gf of 
the UV LF of LBGs at that redshift (i.e. $^Xb = Gf^Hc)- 
We estimate Gf by comparing the model predictions of UV 
LF of LAEs with the observed data. 



As a fiducial model we have chosen n — 1.0 at all red- 
shifts and for Magn, 0.8 x IO^^Mq, 1.5 x IO^^Mq at redshifts 
3 and 4. For higher redshi fts we do not req uire any AGN 
feedback. As in SSS07 and Ijose et al.l (|2013l ). we find that 
with these parameters the observed luminosity functions is 
well reproduced for f*/rj of 0.042, 0.038, 0.037 and 0.044 
respectively for 2: = 3, 4, 6 and 7. We note that the values of 
f*/ri is nearly constant over the redshift range considered. 
These values are tabulated in Table [T] If we use the averag e 
value oiri at z = 3— 7 estimated from lGonzalez et al.l (|2012r ): 
iReddv et all (|2012l '). we get /* = 0.13-0.05 in the same red- 
shift range. This indicates an increase in the baryon fraction 
that is being converted into stars with time. In the top panel 
of Fig. [1] we show our model predictions of best fit UV LF 



of LBGs in dashed lines. The corresponding observational 
data is shown as solid squares. 

The best fit values of Gf that fits the UV LF of LAEs 
are also tabulated in Table [l] Since there is no available 
UV LF data of LAEs at 2 = 6.6 we have fixed Gf to be 
0.5 at that redshift [j . Our model predictions of UV LF of 
LAEs are shown in the top panel (in solid black lines) of 
Fig. [T] along with the observational data. Given Gf one can 
obtain the Lyman-a LF of LAEs by varying f,f csc°' as a free 
parameter. In this way, the observed Lya LFs of LAEs are 
well reproduced for f,f^^r = 0.066,0.040,0.031 and 0.029 
respectively at z — 3.1, 3.8, 5.7 and 6.6. These best fit /./c^c" 
are tabulated in Table [T] We also show in the lower panel of 
Fig[T] our models predictions of the best Lya LFs of LAEs 
along with observed data. Using the /* determined earlier 
we have obtained fasc" ~ 0.5 — 0.6 in the redshift range z — 
3.1 — 6.6. Thus our models predicts no strong evolution in 



^ Note that for predicting UV LF of LAEs at a given redshift we 
use the ft/r] constrained by the observed UV LF of LBGs at a 
slightly different redshift. For example to predict the UV LF of 
LAEs at 2 = 3.1 we use the /*/?? obtained using 2 = 3 observed 
UV LF of LBGs. 
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fcsc°' fro m z = 3.1 to 6.6 wh i ch is consistent wit h the earlier 
studie s of iHaves et al.l(|2010l ): [Blanc et al.l (|201l[ ): IOno etHI 
l|201Cl ). who also found no clear evolution of f^sc" in the 
above redshift range. 



3 TWO POINT CORRELATION FUNCTIONS 
OF LAES 

In order to compute the galaxy-galaxy correlation func- 
tion on all scales one requires the full knowledge of the 
halo occupation distribution (HOD), which describes the 
conditional probability P{N\AI) for A'^ galaxies of a given 
type to reside inside a hal o of mass M (|Bullock et al.ll2002l : 
iBerlind fc Weinberej 120021 ). We modify the phy sically mo- 
tivated HOD models of LBGs (|jose et al.ll2013l . here after 
J13) to compute the same for LAEs that are brighter than 
a threshold Luminosity L^j^" . The mean number of galax- 
ies inside a dark matter halo can be written as Ng{M) — 
fcen{M) + Ns{M) where fcen{M) and Ns{M) are respec- 
tively the mean number of central and satellite galaxies. We 
also assume the central galaxy to be situated at the center 
of the halo and satellite galaxies foll owing NFW density d is- 
tribution dark matter matter halos (jKravtsov et al.ll2004l ). 



3.1 The central galaxy occupation 

In our model the luminosity L of a galaxy initially increases 
and then decreases with its age. Therefore, a galaxy of mass 
M collapsed at two different redshifts can shine with the 
same luminosity L at a given observational redshift. De- 
pending on the redshift of formation these galaxies will be 
either in their increasing phase or decreasing phase of lumi- 
nosity. In between these two epochs (ie two redshifts) the 
galaxy will shine with luminosity greater than L. For such 
a model, the mean occupation number of central galaxies 
with a Lyman-a luminosity above a threshold luminosity 
(or magnitude) Lj^^" inside a halo of mass M at any red- 
shift 2 as in Eq. (15) of J13 is 



J dz, 



Ur.{LT,M,z) = Gs- 



dn{M, Zc) 
dzc 



n{M,z) 



(5) 



Rerezi{M,L 



Lya 



and 2:2 (M, L^^"" , z) are the two redshifts 



at which a galaxy of mass M is to be formed, so that it shines 
with an observed luminosity L^^^" at z. Note a galaxy of 
mass M that has collapsed between these two redshifts z\ 
and 22 will shine with a luminosity greater than L^^"' . We 
show in Fig [2] the halo mass M that can host a galaxy of 
apparent magnitude, rUj^g — 25.3 (of the Lyman-a line) 
at 2 = 3.1. Here, m^g is the narrow band apparent mag- 
nitude and can be computed from L^^" for a fixed set of 
cosmological parameters. The figure clearly shows that ha- 
los of any mass (above a minimum mass Mmin) formed at 
two different redshifts can shine with the same magnitude at 
2 = 3.1. For example, when M — lO^^M© (shown by dotted 
red horizontal line ) these two redshifts are approximately 
3.13 and 5.9 For more details on the evaluation of 21 and 22 
see Section (3) and Fig. (3) of J13. 

We also note that when M < Mmin{L^i^") the central 
galaxy occupation drops to zero, as galaxies below this mass 
scale can never produce a luminosity greater than L^^" . For 




Figure 2. Masses of the dark matter halos that shine with an 
apparent magnitude m^^^ = 25.3 at z = 3.1 as a function of their 
formation (collapse) redshift. Halos of a particular mass formed 
at two different redshifts can shine with the same brightness at 
z = 3.1 which we refer as 21 and 22- For example, when M = 
10^^ Mq (shown by dotted red horizontal line) these two redshifts 
are approximately 3.13 and 5.94. The figure also shows in red 
horizontal solid line the minimum mass of the galaxy that can 
produce a magnitude 25.3 at 2 = 3.1 which is roughly 2.4 X 



771 ^g" = 25.3 and 2 = 3.1 we have shown this threshold mass 
in Fig. [21 We have tabulated Mmi„(L(^") in Table[T]at each 
redshift. From the Table it can be seen that the Mmin is of 
the order of few times lO^^M© at all redshifts and for all 
threshold magnitudes. These threshold masses of galaxies 
that can host an LAEs satisfying the threshold luminosities 
are roughly one order in magnitude smaller than lowest mass 
of LBGs with comparable threshold luminosities in UV light 
(see Table 3 of J13). 

In Fig. [3] we have plotted as thin lines, the average oc- 
cupation number of central galaxies {fcen{M)) as function 
of the mass of the parent halo, calculated using the above 
prescription at different redshifts. For computing the fcen 
we have used G/ and f* fasc" which best fits the UV LF and 
and Lya LF of LAEfl 



3.1.1 The satellite galaxy occupation 

Following J 13 we have used the progenitor mass function 
and our star formation prescription to obtain an estimate of 
the number of subhalos and thereby the number of satellite 
LAEs hos ted by a halo of m ass M. The conditional mass 
function (jMo fc White! 119961 ) gives the comoving number 
density of subhalos of mass Ms which formed at z^ inside a 
region containing a mass M (or comoving volume V) that 



^ Note our LF calculations ignores the presence of satellite galax- 
ies in the parent dark matter halo. However J13 has shown that 
the inclusion of their contribution produces very minor changes 
to the best fitted /»/'/y. 
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have a density contrast 5 ai Zp < Zs. To begin with, J13 
consider regions of mass M which have aheady formed viri- 
ahzed halos by the observational redshift z. In this Umit the 
conditional mass function gives the mass function of subha- 
los in the mass range Ms and Ma + dMs at Za inside a halo of 
mass A4 that collapsed at z < Zp < Zs- The time derivative 
of this conditional mass function is taken as the formation 
rate of subhalos of mass Ma inside a big halo of mass M. 
This is because as in J13 we assume that even though dark 
matter halos are created and destroyed inside the over dense 
volume of mass M, the satellite galaxies formed in these sub- 
halos have survived and can be observed. The star formation 
models for the satellite galaxies are taken to be the same as 
that of central galaxies (given by Eq. ([T])). This is similar to 
the assumption used by previous auth ors (JLee et al.ll2009l : 
iBerrier et a"Lll2006l : IConrov et aLll2006l ') where feedback due 
to halo merging processes are ignored. 

We further assume in our models that no subhalo should 
be formed very close to the formation epoch of parent halo. 
More precisely if t{zp) is the age of the universe when a par- 
ent halo collapsed then all the subhalos formed inside that 
parent halo within a time interval Ato prior to t{zp) do not 
host a satellite galaxy; rather they will be part of the parent 
halo itself. Thus in our models t{zp) — t{zs) ^ Aio where 
t{za) is the age of the universe when the subhalo is formed 
inside the parent halo. The time scale Ato is of the order 
of dynamical time scale of the parent halo. J13 showed that 
this formalism predicts the occupation of satellite galaxies 
Na{M,Lth,z) that can successfully explain the small scale 
angular clustering of LBGs with Ato of the order of dynam- 
ical time scale tdyn of the parent halo. For predicting the oc- 
cupation of Lyman-Q emitting satellites, we again assumed 
Ato ^ tdyn- Interested readers may refer to Section (4.1.2) 
of J13 for more details. 




10" 



10>« 



10" 



10'^ 



Figure 3. The halo occupation distribution, {Ng{M) = 
fcen(M) + Na(M)), as a function of the mass of the hosting 
halo, as predicted by our models at various redshifts and thresh- 
old magnitudes. In each panel the thin curves corresponds to 
fcen(M) where as the thick curves give the total occupation 
Ng(M) = fcen{M) + Ns{M). The solid and dotted lines cor- 
respond to two different threshold magnitudes. The curves are 
obtained using a fiducial set of model parameters (given in Sec- 
tion 2) that reproduce the observed luminosity function. In addi- 
tion to that, to obtain the satellite occupation {Na{M)), we have 
adopted Aio = 2.0tdyn (see text for details). 



3.1.2 The total halo occupation 

We have plotted in Fig. O the total occupation number of 
galaxies, Ng{M) = fcen{M) + Na{M) as a function of the 
mass of the parent halos at various redshifts and for different 
threshold magnitudes (ie ti^^"). For each z we have used 
the best fitted parameters given in Table [T] These curves are 
shown in thick lines and thin lines correspond to occupation 
number of central galaxies. To obtain Na{M), apart from 
the fiducial model parameters that fit the observed luminos- 
ity function, we have adopted Aio ~ "itdyn. These values of 
Ato are chosen as they are later used in Section 3] for ex- 
plaining the small angular scale clustering. We find that for 
each limiting magnitude the mean numbers of central and 
satellite galaxies are monotonically increasing with the mass 
of the hosting halo mass. 

It is of interest to calculate the average value of Ng = 
fcen + Na, which is tlic average number of LAEs occupied 
in dark matter halos above the threshold mass Mmin . As in 
J13 we define this quantity as 

J^ . dM{fceniM,z) + Na{M,z))n{M,z) 



{Ny}{L^^ 



Lya 



J^ dMn{M,z) 



(6) 
We give the value of (Ng) in Table [T] for different magnitude 
thresholds and redshifts. From this table it is clear that at 
2 = 3.1 there is about 1-2 % of LAEs in dark matter ha- 



los above the threshold mass Mmin. This number is much 
smaller than that of LBGs (see Table 3 of J13) which was 
roughly 40%. However as we go to higher redshift, the aver- 
age Ng of LAEs increases. For z ^ 5.7 this is roughly 20 % 
which is roughly 50 % of the mean occupation of LBGs at 
higher redshifts. The mean Ng desc ribed above can be in- 
voked as the duty cycle of LAEs fsee lLee et al.l ((20091), J13) 
hence this directly implies that LAEs have a smaller duty 
cycle compared to LBGs. 

Similarly in Table [l] we also give the average mass Mav 
of a halo hosting the LAEs, for each redshift and luminosity 
threshold. This average mass is defined as. 



Mav — 



Jp°° dMM {fcenjM, Z) + NajM, z)) UJM, z) 

j;^dMUn{M,z)n{M,z) 



(7) 



One can find from Table [T] that the average mass of LAEs 
of the order of 10^° - lO^M© at all redshifts and threshold 
magnitudes we are considering. For example, at z = 3.1 and 
for m^^° = 25.3 we find that Mav = 7 x 10^" M© and at 
2 = 5.7 for m^^°^ = 25.5 Mav = 1.4 x IO^Mq. We also 
note that, at a given redshift Mav of LAEs increases with 
higher threshold luminosity. Further, the average masses 
halos, hosting LBGs of similar continuum luminosity, are 
higher by an order of magnitude (given in Table 3 of J 13). 
This is in agreement with the suggestions that LAEs are re- 
sidin g in smaller mass halos compared to LBGs (|Quchi et al.l 
l2010t ). 
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z 


f*h 


Gf 


/. pLya 
J*Jesc 


Lya 


log(LLy^) 


M„j„/Mq 


Mav/MQ 


N, 


Ug 


h 


bgiO) 


~ 3.1 


0.042 


0.041 


0.066 


25.3 


42.1 


2.4 X 10^0 


7.1 X 10^0 


0.011 


1.8 


2.1 


1.5 ± 0.7 










24.8 


42.3 


3.8 X 10^° 
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0.014 


1.1 
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0.059 
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24.7 
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10" 


2.2 X 10^1 


0.024 


0.4 


3.5 


2.8 ± 0.5 
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0.037 
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26.0 


42.4 


5.0 X lO^o 


9.1 X lOio 


0.190 


1.8 


5.1 


5.5 ± 0.4 










25.5 


42.6 


7.9 X lO^o 


1.4 X 10" 


0.180 


0.8 


5.6 


6.1 ± 0.7 


~ 6.6 


0.044 


0.50 


0.029 


26.0 


42.4 


4.3 X 10^0 


7.5 X 10"' 


0.230 


1.1 


6.1 


3.6 ± 0.7 










25.5 


42.6 


6.8 X 10^" 


1.2 X 10" 


0.213 


0.5 


6.7 


6.0 ± 2.2 



Table 1. Best fitted and derived parameters of our models. (1) Redshift; (2) ft/r], the parameter related to the hght to mass ratio of 
galaxies; (3) Gf, the fraction of LBGs that shows detectable Lyo emission; (4) ftfeJi"', the parameter related to the net Ly« emission 
from LAEs; (5) the limiting magnitude of LAE samples; (6) approximated limiting Lya luminosity corresponding to m^^; (7) minimum 
mass of dark matter matter halos that can host LAEs; (8) average mass of dark matter matter halos that can host LAEs; (9) average 
occupation of LAEs in dark matter halos; (10) number density of LAEs in units o f 10~ (h/Mpc) ; (11) large scale galaxy bias of LAEs 
of the given threshold magnitude; (12) best estimate of bias bv lOuchi et alj 1120101 ). 



3.2 The correlation functions. 

In the framework of ha lo model, the total cor relation func- 
tion can be written as (jCoorav fc Shethll2002l ) 



UR) = aR,^)+C(R^^) 



(8) 



where each term on RHS has contributions from central as 
well as satellite galaxies. On scales much bigger than the 
virial radius of a typical halo, the clustering amplitude is 
dominated by correlations between galaxies inside separate 
halos. On the other hand, on scales smaller than the typical 
virial radius of a dark matter halo, the major contribution 
to galaxy clustering is from galaxies residing in the same 
halo. 

The 2 halo term, ^^''(''i ^y Lth'^)y ^^ galaxies with lumi- 



nosity greater than L^^ 
(|Peeblej|l980l ). 



Cs iR,z,L^ 



Lya 



at z, can now be calculated using 



dfc 2Sin(fc^ 2h/, r Lya\ 

2^^" kR ^« ^'='^'^*h ) 



(9) 



Here Pg^(k, L^^" , z) is the luminosity dependent galaxy 



power spectrum, which is computed as 

Pg^{k,Lf^",z) = bl{k,Lf^",z)Pii„{k,z). 



(10) 



In above equation Piin{k, z) is the linear power spectrum of 
CDM density fluctuations. The luminosity dependent galaxy 
bias obtained by adding the contributions of both the central 
and satellite galaxies is, 



bg{k,z,Lf^°') = 



1 



^liLj 



Lya 



\jcen(-lJth 1 



A/, z) 



dMn{M,z)b{M,z) 

-{N4Lfr,M,z)) 
u{k, M, z) 



(11) 



As before, n[M, z) is the ST halo mass function, and 

nl^L'^r. z) = j dMn{M, z)iU„{M, z) + {N4Lfr,M, z)) 

(12) 
is the total number density of galaxies, including both cen- 
tral and satellite galaxies. Further, h{M, z) is the mass de- 
pendent halo bias f actor provided by the fitting function of 
Sheth and Tormen (|Sheth fc Tormenlll999l : ICoorav fc ShethI 
|2002| . also see Eq. (7) of J13). In Eq. ^ u{k,M) is the 



Fourier transform of dark matter density profile normal- 
ized by its mass, i.e u{k,M,z) = p{k,M,z)/M. For the 
present calculations we use the Navarro Frenk White (NEW) 
form for the dark matter density distribution (jNavarro et al.l 
[l997|) . On scales much larger than the virial radius of typi- 
cal halos u{k, M, z) — 1 and hence the luminosity dependent 
galax;y bias given by Eg. 1111 becomes a constant. 

We give in Table [T] this constant galaxy bias computed 
using our models. We note that the large scale galaxy bias 
predicted by our models c ompares well with obser vat ionally 
derived galaxy bias from lOuchi et al.l (|2010|) (given in the 
last column of Table [l]) at each redshift and for nearly all 
threshold luminosities. We also find that the galaxy bias 
for a given z increases with increasing luminosity and is 
systematically higher for high redshift galaxies. Further, the 
large scale bias, computed here for LAEs of a given Lyman 
a threshold luminosity, is smaller than that LBGs with the 
same continuum luminosity (see Table 2 of J 13 ) . This result 
has been observed previously by lOuchi et al.l (|201Cl ) . This 
is because, as we have seen in the last section, LAEs of 
given Lyman a luminosity are residing in dark matter halos 
of smaller mass compared to LBGs of the same continuum 
UV luminosity. This directly results in smaller clustering 
strength of LAEs compared to LBGs. We elaborate on this 
interesting fact in the next section. 



3.3 The total correlation functions 

The 1-halo term ^l'^(R) uses the standard assumption that 
radial distribution of satellite galaxies inside a parent halo 
follow the dark matter density distribution (NFW profile) 
(ICoorav fc Shethll2002l) . In this case the 1-halo te rm is given 
by (|Tinker et all (|2005l ): ICoorav fc Ouchil (|2006l )) 



e'iLfr, R, z) = ^ / dM n{M, z) x [2/_(Lf,-", M, z) 

{N.\Lfr,M,z) P-^^^f\^'K 

(iV.(iV.-l)|ifr,A/,->^^^^^f^].(13) 

Here, as before, n'^{L^^",z) is the total number density 
of galaxies which includes both the central and satellite 
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Figure 4. The angular correlation function of LAEs, taking into account the one halo and the two halo terms, at redshifts 3.1, 3.8, 5.7 
and 6.6 for various limiting magnitudes. Each row corresponds to a particular redshift, which is labelled in the first panel of that row. 
In each panel clustering predictions for galaxies with a given threshold magnitude is shown. Our fiducial model predictions of galaxy 
angular correlation functions with Ato = "itdyn is shown in solid black lines. The dashed red curves arc our model predictions assuming 
a sub-Poisson distributio n for galaxies insid e dark matter halos (see text for details). The data points and error bars shown by black 
open rectangles are from lOuchi et al.l l|201Clf ). 



galaxies. F urther pnfw is the NFW profile of dark mat- 
ter density (JNavarro et al.lll997l ) inside a collapsed halo and 
AjvFvi/ (r|M) is the conv olution of this density profile with 
itself IjSheth et al.l 1200 If ). For the halo conce ntration pa- 
ramet er we use the fitting functions given by IPrada et al.l 
(|2012l ) (Eq. (14-23) of their paper). Also to begin with, fol- 
lowin g the N-body simulat i ons and sem i -anal ytical models 
(eg. iKravtsov et al.l 120041 : IZheng et al.ll2005h , we assume 
that the number of satellites inside a parent halo forms a 
Poisson distribution. Thus we have {Ns{Ns — 1)) ~ N"^ , 
when Ng is large. 



Finally we compute luminosity dependent angular cor- 
relation function w{6, z) from the spatial correlation func- 



tion using Limber equation ijPeebleslllQSd . J 13) 



w{e,z) = / dz' N{z) / dz" N(z")£_g {z,r{6\z ,z")) 
Jo Jo 

(14) 
where r{9;z',z") is the comoving separation between two 
points at z' and z" subtending an angle 9 with respect to 
an observer today. Here we have also incorporated the nor- 
malized redshift selection fun ction, N(z), of the obser ved 
population of galaxies given in lOuchi et al.l (|2010l . 120081 ). In 
Eq. (|14|l we neglect the redshift evolution of clustering of the 
galaxies detected around z. Hence the spatial two-point cor- 
relation function ^g (r, z) is always evaluated at the observed 
redshift. 
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agn and K are kept to their fiducial 

AU other 



4 COMPARISONS WITH OBSERVATIONS 

The total angular correlation functions computed using our 
prescription for four redshifts and two threshold magnitudes 
are overplotted on the observed data in Fig. |4l The cosmo- 
logical parameters, M, 

value. AU these models assumes Aio = 2 tdyn 
other model parameters are obtained by fitting the LFs and 
are given in Table [T] Firstly, we note that angular correla- 
tion functions predicted by our model (for 6 ^ 10") is in very 
good agreement with the observed data at all redshifts and 
for all threshold luminosities. The amplitude of these angu- 
lar correlation functions are higher than that of LBGs (see 
Fig. 9 of J 13) mainly because LAEs have a much narrower 
redshift distribution function N{z). 

We note that our models systematically overpredicts 
clustering at very small angular scales {6 SC 10") at all red- 
shifts and for all magnitudes. That is the contribution of one 
halo term predicted by our models to the w(S) for small val- 
ues of 9 is not distinctly evident from the observed w{6). We 
now investigate whether this discrepancy could be due to our 
earlier assumption that galaxies are Poisso n distributed in 
dark m atter halos. N-body simulations (eg. iKravtsov et al.l 
(I2OOJ)) showed that the HOD follows a Poisson distribu- 
tion when the total number of gala xies inside a halo exceeds 
unity. Number of earlier w orks ( Scocci marro et al.l I2OOII : 
IWechsler et al.l 1200 ll : ISullock ct al. 2002) have pointed out 
that, galaxies follow sub-Poisson distribution when Ng ^ 1. 
Note from FigO the average halo occupation {{Ng}) of LAEs 
(3.1 ^ z ^ 6.6) in dark matter halos are ty pically smaller 
than unity. The refore we choose the model of lBuUock et al.l 
(|2002l ) (see also lHamana et al.l (|200J)) for the average num- 
ber of galaxy pairs inside a dark matter halo of mass M, 



{Ng{Ng-l)) = 



N^ ln{ANg)/ln{4) 




0.25 < iVj < 1 (15) 



Ng ^ 0.25 



In this case the one halo term is computed as in 
iHamana et al.l (|2004l '). where 



e\L 



Lya 
th 



,R,z)^ 



k: 



{Ng{Ng-l)\Lfr,M,z) 



dM n{M, z) 
\nfw{R\M) 



(16) 



In Fig. U we show our model predictions (in red dashed 
lines) when one uses the above sub-Poisson distribution for 
galaxy occupation. One can see from the figure that our 
model predictions of angular correlation functions are con- 
sistent with the observed data even at very small angular 
scales {6 ^ 10") for all redshifts and threshold magnitudes. 

As discussed earlier, galaxies selected by their thresh- 
old AB magnitude in Lyman a luminosity, are less strongly 
clustered than the galaxies selected by same UV threshold 
AB magnitude (See Fig. 9 and Table 2 of J13 also). This 
can be understood in the following way. In Fig. [5] we have 
plotted the difference in the observed Lyman a and UV mag- 
nitudes of LAEs at 2 = 3.1 against th eir Lyman a magni- 
tudes. The data points are from from lOuchi et al.l (|2008n . 
The figure clearly shows that the observed Lyman a appar- 
ent magnitude of LAEs residing in a halo of given mass are 
higher than their continuum magnitude roughly by a factor 



Si 2 



- 




Lya-. 



of AruAB = rriAB 



Lya 



2. Or equivalently, the Lyman 



Figure 5. A comparison of the Lyman a AB magnitudes {m^g ) 
and UV AB magnitudes {m^^) of observed Lyman a emitters at 
z = 3.1. The Y-axis shows the difference between UV and L yman 
a magnitudes. The data points are from lOuchi et al.l l|2008 



a line luminosity of LAEs are typically higher than their 
continuum UV luminosity. This inturn implies that, LAEs 
of a given Lyman a luminosity, can be found in low mass 
dark matter halos compared to an LBG of the same UV 
threshold magnitude. As a result, the galaxy bias of LAEs 
of a given Lyman a threshold magnitude are smaller than 
LBGs of the same threshold magnitude in UV. 

The above discussions suggest that at z = 3.1, the 
galaxy bias of LAEs with a line threshold magnitude, JTi^g", 
must be comparable to the galaxy bias of LBGs with contin- 
uum threshold magnitude, rnJ^^, where m}^^ ~ ^^ab' + 2- 
In order to demonstrate this more elaborately, we computed 
the galaxy bias of LBGs (using the prescription of J13) at 
2; = 3.1 for threshold magnitude m^]^ = 26.8 and compared 
it with that of LAEs of threshold magnitude rrij^'^ = 24.8 
(here Atuab = 2). Interestingly we found that at 2; = 3.1 
galaxy bias of LBG samples (using the best fitted parameters 
of LBG LF) with tti'ab = 26.8 is bg = 2.35 which compares 
very well with galaxy bias oi bg = 2.2 that we obtain for 
LAEs with m^g = 24.8. We performed similar exercise at 
other redshifts and found comparable galaxy bias for LAEs 
and LBGs when appropriate limiting magnitude are chosen. 
This implies that the connection between dark matter halo 
and UV luminosity found by J13 for the LBGs is valid for 
the LAEs as well. The lower clustering amplitude one finds 
for LAE samples can be attributed to the narrow band tech- 
nique having higher efficiency in picking up low luminosity 
LBGs compared to the LBG surveys. 



5 DISCUSSION AND CONCLUSIONS 

We have presented here a physically motivated semi- 
analytical model of galaxies to understand the clustering 
of high redshift Lyman-a emitters. In our models we as- 
sign luminosities (UV and Lya) to galaxies residing in dark 
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matter halos by a physical model of star formation. The 
free parameters of this star formation model are then con- 
strained by the observed high redshift luminosity functions. 
Further using the semi analytical approach of J13, we com- 
bined these constrained models of star formation with the 
halo model and predicted the angular correlation functions, 
w{9) of LAEs. Our model predictions of w{0) of LAEs with 
different threshold magnit udes compares rem arkably well 
with the observed data of lOuchi etall (|2010l ') for 9 > 10" 
in the redshift range z ~ 3 — 7. 

Our models assumed that only a fraction G/ of halos 
hosting LBGs host LAEs. The low amplitudes of LPs LAEs 
compared to that of LBGs at various redshifts suggest that 
Gf is typically smaller than unity. We found that Gf ranges 
from 4% — 50% in the redshift range 3.1 — 6.6. Such small 
values of Gf directly result in a low halo occupation for 
LAEs compared to LBGs. Our models for LBG halo occu- 
pation (see J 13) suggested that about 40% dark matter halos 
above a minimum mass Mmin host LBGs above a threshold 
luminosity. On the other hand, the average number LAEs 
occupying dark matter halos, {Ng), are much smaller, rang- 
ing from 1% — 23% in the above redshift range. As low halo 
occupation is related to duty cycle (see J13), we can con- 
clude that Lyman-a emission has much smaller duty cycle 
compared to the UV emitting phase of the high redshift 
galaxies. 

We find that, for observationally d etermined r] from 
iBouwens et al.l (|2012l ): lReddv et al.l (|2012h studies at various 
redshifts, the escape fraction of Lyman-a from LAEs varies 
from ~ 0.5 to ~ 0.6 for the redshift range 3.1 ^ z ^ 6.6. 
Thus in our studies /i^c" is not strongly evolving in the 
above redshift range and such a trend is consistent with 
previous studies. The lack of sharp decrease in /esc" in this 
redshift range is consistent with the epoch of reionization 

Zre ^ 6.6. 

We found that the average masses, Mav of halos hosting 
LAEs of a given threshold Lyman-a luminosity are smaller 
than that of LBGs with same continuum threshold luminos- 
ity. For Lyman a apparent magnitudes 25-26, Mav ranges 
typically from 7.0 x 1O^°M0 - 1.2 x 10^^ Af© and are smaller 
than typical LBGs of similar continuum magnitude by a fac- 
tor 10. Because of their smaller halo masses LAEs of a given 
threshold line luminosity, have smaller galaxy bias and are 
more weakly clustered than LBGs of same continuum lumi- 
nosities. Our model predictions of galaxy bias c ompares very 
well w ith observationally derived galaxy bias of lOuchi et al.l 
l|2010f ). We also found that, the smaller masses and biases 
of LAEs of given line luminosity compared to LBGs of sim- 
ilar continuum luminosity is the simple consequence of the 
fact that line luminosity of LAEs are typically higher than 
their continuum UV luminosity. In particular. Fig. [S] shows 
that the difference between the observed continuum magni- 
tude, m^AB^ a-nd line magnitude, rrij^g , of LAEs at 2 = 3.1 
is AmAB ~ 2. We further showed that the galaxy bias of 
LBG sample with m^AB = 26.8 is almost same as the galaxy 
bias of LAEs with m^g — 24.8. This clearly suggests that 
LAEs belong to the same galaxy population of LBGs with 
narrow band technique having more efficiency in picking up 
galaxies with low UV luminosity or low mass. 

Finally we note that our simple model, that uses Pois- 
son distribution for satellite galaxies in dark matter halos, 
predicts an excess in the small angular scale [Q ^ 10") clus- 



tering. The angular correlation functions on these scales are 
dominated by one halo term. We demonstrate that this dis- 
crepancy can be removed if one uses a sub-Poisson distri- 
bution for galaxies inside dark matter halos when the total 
galaxy occupation Ng is smaller than unity. Since G/ is al- 
ways smaller than unity, our models predict a low occupation 
of LAEs in dark matter halos (See Fig[3]) for all redshifts and 
for various threshold magnitudes. The total occupation, Ng 
is bigger than unity only for very massive dark matter halos. 
The small angular scale clustering of LAEs can in principle 
be used constrain the statistics of galaxy occupation in dark 
matter halos when the mean occupation, Ng, is smaller than 
unity. 

We have showed tentatively that the clustering of LAEs 
can be understood in the framework of halo model and phys- 
ical models of galaxy formation. A clear understanding of 
these LAEs requires the complete knowledge of Lyman-a 
radiative transport in the interstellar medium. Nevertheless, 
our simple model suggests that, LAEs belong to the same 
galaxy population as LBGs and narrow band techniques are 
picking up galaxies with strong Lyman-a lines. 
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